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Abstract The purpose of this paper is to provide a mathematical analysis of the 
Adler- Wiser formula relating the macroscopic relative permittivity tensor to the mi- 
croscopic structure of the crystal at the atomic level. The technical level of the pre- 
sentation is kept at its minimum to emphasize the mathematical structure of the 
results. We also briefly review some models describing the electronic structure of fi- 
nite systems, focusing on density operator based formulations, as well as the Hartree 
model for perfect crystals or crystals with a defect. 



1 Introduction 

Insulating crystals are dielectric media. When an external electric field is applied, 
such an insulating material polarizes, and this induced polarization in turn affects 
the electric field. At the macroscopic level and in the time-independent setting, this 
phenomenon is modelled by the constitutive law 

D = e e M E (1) 

specifying the relation between the macroscopic displacement field D and the 
macroscopic electric field E. The constant £o is the dielectric permittivity of the 
vacuum, and £m the macroscopic relative permittivity of the crystal, a 3 x 3 sym- 
metric tensor such that £m > 1 in the sense of symmetric matrices (k^ejyjk > |k| 2 for 
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all k G R 3 ). This tensor is proportional to the identity matrix for isotropic crystals. 
Recall that D is related to the so-called free charge pf by the Gauss law div(D) = pf 
and that the macroscopic electric field E is related to the macroscopic potential V 
by E = — VV, yielding the macroscopic Poisson equation 

-div(e M VV)=p f /£b. (2) 
In the time-dependent setting, (HJ becomes a time-convolution product: 

D{r,t)=e Q J Eu{t-t')E{r,t')dt'. (3) 

Fourier transforming in time, we obtain 

TD(r,(o) = >FE M {(o)TE(r,(o), 

where, as usual in Physics, we have used the following normalization convention 
for the Fourier transform with respect to the time-variable: 

y/(r,a>) = J + f(r,t)e im 'dt 

(note that there is no minus sign in the phase factor). The time-dependent tensor 
£m in 03 can be seen as the zero-frequency limit of the frequency-dependent tensor 

Of course, the constitutive laws fl} (time-independent case) and (0) (time- 
dependent case) are only valid in the linear response regime. When strong dielectric 
field are applied, the response can be strongly nonlinear. 

The purpose of this paper is to provide a mathematical analysis of the Adler- 
Wiser formula JT] [36) relating the macroscopic relative permittivity tensor £m (as 
well as the frequency-dependent tensor f e M (a))) to the microscopic structure of the 
crystal at the atomic level. 

In Section[2] we discuss the modelling of the electronic structure of finite molec- 
ular systems. We introduce in particular the Hartree model (also called reduced 
Hartree-Fock model in the mathematical literature), which is the basis for our anal- 
ysis of the electronic structure of crystals. This model is an approximation of the 
electronic iV-body Schrodinger equation allowing to compute the ground state elec- 
tronic density of a molecular system containing M nuclei considered as classical 
particles (Born-Oppenheimer approximation) and N quantum electrons, subjected 
to Coulomb interactions. The only empirical parameters in this model are a few 
fundamental constants of Physics (the reduced Planck constant h, the mass of the 
electron m e , the elementary charge e, and the dielectric permittivity of the vacuum 
£o) and the masses and charges of the nuclei. In this respect, this is an ab initio, or 
first-principle, model in the sense that it does not contain any empirical parameter 
specific to the molecular system under consideration. 
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We then show, in Section [3] how to extend the Hartree model for molecular sys- 
tems (finite number of particles) to crystals (infinite number of particles). We first 
deal with perfect crystals (Section [3. 2\ , then with crystals with local defects (Sec- 
tion l3.31 l. The mathematical theory of the electronic structure of crystals with local 
defects presented here (and originally published in Q) has been strongly inspired 
by previous works on the mathematical foundations of quantum electrodynamics 
(QED) [18, 20, 19 1 . In some sense, a defect embedded in a insulating or semi- 
conducting crystal behaves similarly as a nucleus embedded in the polarizable vac- 
uum of QED. 

In Section |4] we study the dielectric response of a crystal. First, we focus on 
the response to an effective time-independent potential V, and expand it in powers 
of V (Section [4. U . The linear response term allows us to define the (microscopic) 
dielectric operator e and its inverse e _1 , the (microscopic) dielectric permittivity 
operator, and also to define a notion of renormalized charge for defects in crystals 
(Section [4.21 . In Section 1431 we derive the Adler- Wiser formula from the Hartree 
model, by means of homogenization arguments. Loosely speaking, a defect in a 
crystal generates an external field and thereby a dielectric response of the crystal. 
If a given local defect is properly rescaled, it produces a macroscopic charge (cor- 
responding to the free charge pf in (O and the total Coulomb potential converges 
to the macroscopic potential V solution to (f2]i where £m is the tensor provided by 
the Adler- Wiser formula. A similar strategy can be used to obtain the frequency- 
dependent tensor f £m(©) (Section l4~4l . 

As trace-class and Hilbert-Schmidt operators play a central role in the mathemat- 
ical theory of electronic structure, their definitions and some of their basic properties 
are recalled in Appendix for the reader's convenience. 

The mathematical results contained in this proceeding have been published Q[8] 
9 1, or will be published very soon IflOll . The proofs are omitted. A pedagogical effort 
has been made to present this difficult material to non-specialists. 

As usual in first-principle modelling, we adopt the system of atomic units, ob- 
tained by setting 

h = 1 , m e = l, e = 1 , = 1 , 

4ne 

so that reads in this new system of units: 

-div(e M VV) =4%pf. (4) 

For simplicity, we omit the spin variable, but taking the spin into account does not 
add any difficulty. It simply makes the mathematical formalism a little heavier. 
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2 Electronic structure models for finite systems 



Let 3^ be a Hilbert space and (-| ) its inner product (bra-ket Dirac's notation). Recall 
that if A is a self-adjoint operator on ffl and <j> and \j/ are in D(A), the domain 
of A, then (0|A|y/) := (0|Ayf) = (A(j>\y). If A is bounded from below, the bilinear 
form (<p,y) i-> (0 |A | y/) can be extended in a unique way to the form domain of A. 
For instance, the operator A = —A with domain D(A) = H 2 (M. d ) is self-adjoint on 
L 2 (R d ). Its form domain is H l (R d ) and (<j)\A\ yr) — j Rd V0 ■ Vy/. In the sequel, we 
denote by ^(Jff) the vector space of bounded self-adjoint operators on Jtf?. 

For k = 0, 1 and 2, and with the convention H°(M?) = L 2 (R 3 ), we denote by 

N 

/\H k (R 3 ) := { V € | f(r Kl)! - • • ,r p(w) ) = e(p)«P(n,- • • Vp e §n} 

(=i 

(where is the group of the permutations of {1, • • • ,A^} and e{p) the parity of p) 
the antisymmetrized tensor product of spaces H k (M. 3 ). These spaces are used to 
describe the electronic state of an electron system. The antisymmetric constraint 
originates from the fact that electrons are fermions. 



2.1 The N-body Schrddinger model 

Consider a molecular system with M nuclei of charges z\, — - ,Zm- As we work in 
atomic units, Zk is a positive integer. Within the Born-Oppenheimer approximation, 
the nuclei are modelled as classical point-like particles. This approximation results 
from a combination of an adiabatic limit (the small parameter being the square root 
of the ratio between the mass of the electron and the mass of the lightest nucleus 
present in the system), and a semi-classical limit. We refer to JUG] ar, d references 
therein for the mathematical aspects. 

Usually, nuclei are represented by point-like particles. If the M nuclei are located 
at points Ri , • • • , Rm of M 3 , the nuclear charge distribution is modelled by 

M 

k=l 

where 8r k is the Dirac measure at point R^. The Coulomb potential generated by 
the nuclei and seen by the electrons then reads 

M 

y nuc (r) = y — - — 
U fc|r-R*| 

(the minus sign comes from the fact that the interaction between nuclei and electrons 
is attractive). In order to avoid some technical difficulties due to the singularity of 
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the potential generated by point-like nuclei, the latter are sometimes replaced with 
smeared nuclei: 

M 

P nuc (r)-Eza(r-R^), 

k=l 

where % is a smooth approximation of the Dirac measure 8q, or more precisely a 
non-negative smooth radial function such that L 3 % = 1, supported in a small ball 
centered at 0. In this case, 

y „u C(r) := _ (p „uc^| . r i )(r) = _r ppil dr > 

is a smooth function. We will sometimes denote this smooth function by Vpnuc in 
order to emphasize that the potential is generated by a non-singular charge distribu- 
tion. 

The main quantity of interest in our study is the electrostatic potential gener- 
ated by the total charge, which is by definition the sum of nuclear charge p nuc and 
the electronic charge p el . According to the Born-Oppenheimer approximation, elec- 
trons are in their ground state, and p el is a density associated with the ground state 
wavefunction %. Let us make this definition more precise. 

Any (pure) state of a system of N electrons is entirely described by a wavefunc- 
tion f G /\^ =1 L 2 (R 3 ) satisfying the normalization condition ||f H^^at) = 1. The 
density associated with f is the function p^ defined by 



Clearly, 



p^(r) = N / |f(r,r 2 ,--- ,r N )\ 2 dr 2 - ■ -dr N . 

Pf>0, p^eL'fM 3 ), and / py=N. 

Js? 



(5) 



It can be checked that if f e A£i^'(K 3 ), then ^p g i/^R 3 ), which implies in 
particular that p w € L l (R 3 ) n L 3 (M 3 ) . 

The ground state wavefunction % is the lowest energy, normalized eigenfunction 
of the time-independent Schrodinger equation 



N 



h n v = ev, f e/\i/ 2 (R 3 ), ||f|| i2(RW) = i, (6) 

where Hn is the electronic Hamiltonian. The latter operator is self-adjoint on 
/\f =l L 2 {R 3 ), with domain /\? =l H 2 (R 3 ) and form domain /\f =l H l {R 3 ), and is de- 
fined as 

^ = -^E^ + Lv nuc (r,)+ I rr 1 ^- (7) 

i !=l i=l l<i<j<N' ' l\ 

The first term in the right-hand side of models the kinetic energy of the electrons, 
the second term the Coulomb interaction between nuclei and electrons and the third 



6 



Eric Cances, Mathieu Lewin and Gabriel Stoltz 



term the Coulomb interaction between electrons. For later purposes, we write 

H N = T + V n e.+ Ve«, 



where 



i£>„ Ke = I> nttc (r«), V ( 



i=l l<i<j<N l r ' 



It is proved in [37| that if the molecular system is neutral (Y!k=i^k = N) or posi- 
tively charged (Y,f=i Zk > N), then the essential spectrum of Hn is an interval of the 
form [Lpj,+°°) with En < and < if N > 2, and its discrete spectrum is an 
increasing infinite sequence of negative eigenvalues converging to En- This guaran- 
tees the existence of ffo. If Eq, the lowest eigenvalue of Hn is non-degenerate, % 
is unique up to a global phase, and p el = p^o is therefore uniquely defined by (f5]). 
If Eq is degenerate, then the ground state electronic density is not unique. As the 
usual Born-Oppenheimer approximation is no longer valid when Eq is degenerate, 
we will assume from now on that Eq is a simple eigenvalue. 

Note that % can also be defined variationally: It is the minimizer of 



;V 



MU¥\H N \¥), f e/\H l (R 3 ), \\nm&») = lj- (8) 

Otherwise stated, it is obtained by minimizing the energy (f |//Ar|f) over the set of 
all normalized, antisymmetric wavefunctions f of finite energy. 

Let us mention that, as in the absence of magnetic field, the A^-body Hamiltonian 
is real (in the sense that it transforms a real-valued function into a real-valued func- 
tion), there is no loss of generality in working in the space of real-valued iV-body 
wavefunctions. Under the assumption that Eq is non-degenerate, ^ has exactly two 
minimizers, % and — ffo, both of them giving rise to the same electronic density. 



2.2 The N-body Schrodinger model for non-interacting electrons 

Neither the Schrodinger equation (O nor the minimization (JHJ can be solved with 
standard numerical techniques when exceeds two or three. On the other hand, 
these problems become pretty simple when the interaction between electrons is ne- 
glected. In this case, the A^-body Hamiltonian is separable and reads 

N 1 

H N = T + V m = Y,K where h Ti = --A T . +V nuc 

i=i l 

is a self-adjoint operator on L 2 (M 3 ) with domain H 2 (M?) and form domain H 1 (M 3 ), 
acting on functions of the variable r,. It is known that the essential spectrum of 
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h is [0, +°o) and that the discrete spectrum of h is an increasing infinite sequence 
of negative eigenvalues converging to 0. Let us denote by £1 < £2 < £3 < • • ■ the 
eigenvalues of h counted with their multiplicities (it can be shown that £1 is simple) 
and let (0,-);>o be an orthonormal family of associated eigenvectors: 



h<j>i = Ei(j>i, £1 < £2 < £3 < 



0e// 2 (K 3 ), 



The eigenfunctions 0; are called (molecular) orbitals and the eigenvalues £, are 
called (one-particle) energy levels. 

It is easy to check that if £ w < £#+1, then 



inf j(f|i/°|f), f e/\H l 



L 2 



= 1 



(9) 



has a unique solution (up to a global phase) given by the Slater determinant 



%(r u --- ,r N ) = 



01 Or) 0i (r 2 ) 
02 (ri) 02(r2) 



0v( r i) 0iv(r 2 ) 



0i(riv) 
02 (rjy) 



0w(rw) 



and that the ground state electronic density (O takes the simple form 



(10) 



p el (r) = ll0<(r 



The above description of the electronic states of a set of /V non-interacting elec- 
trons in terms of orbitals cannot be easily extended to infinite systems such as crys- 
tals (the number of orbitals becoming infinite). For this reason, we introduce a new 
formulation based on the concept of one-particle density operator, here abbreviated 
as density operator. 



2.3 Density operators 

The (one-particle) density operator of a system of /V electrons is an element of the 
convex set 

@ N = {ye .Y(L 2 (R 3 )) |0<7< 1, Tr(y)=/V}. 

Recall that if A and B are two bounded self-adjoint operators on a Hilbert space 3^, 
the notation A < B means that (y|A|yA) < (y\B\y) for all y/ e 

Any density operator 7 € &n is trace-class, hence compact (the basic properties 
of trace-class operators are recalled in the Appendix). It can therefore be diagonal- 
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ized in an orthonormal basis: 

y=X>l<fcX<fcl with (M,) = ay. (id 
1=1 

The eigenvalues n, are called occupation numbers; the eigenfunctions 0, are called 
natural orbitals. The conditions < y < 1 and Tr(y) = N are respectively equivalent 
to 

< Hi < 1 and £ m = N. 

i=i 

The fact that < n, < 1 is a mathematical translation of the Pauli exclusion principle, 
stipulating that each quantum state |0 ; ) is occupied by at most one electron. The sum 
of the occupation numbers is equal to N, the number of electrons in the system. The 
density associated with y is defined by 

p y (r) = £» i |0,(r)| 2 , (12) 

;=i 

this definition being independent of the choice of the orthonormal basis (0,),>i 
in (flTT l and satisfies 

Pr>0, p y g l} (M 3 ), and / p Y = N. 
The kinetic energy of the density operator y is defined as 

ny):=^Tr(|V| 7 |V|), 

and can be finite or infinite. Recall that the operator | V| is the unbounded self-adjoint 
operator on L 2 (M. 3 ) with domain H 1 (M 3 ) defined by 

V0 G H\R 3 ), (^(|V|0))(k) = |k|(JT(0))(k) 

where & is the unitary Fourier transform 

The kinetic energy of a density operator y decomposed as (fTTT i is finite if and only 
if each 0, is ini/'(R 3 ) and E|t!Jn,-|| V0j-|| 2 2( , K 3j < oo, in which case 

7 , (7) = |l"i||V0i||i 2(R3) . 
z i=i 

As |V| is the square root of —A (i.e. |V| is self-adjoint, positive and |V| 2 = —A), 
the element Tr(|V|y|V|) of R+ U {+°°} is often denoted by Tr(-4y). Using this 
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notation, we can define the convex set &n of the density operators of finite energy 

as 

&n = |yGJ^(L 2 (R 3 )) I 0< y< 1, Tr(y) = N,Tr{-Ay) < 00} . 

Lastly, it is sometimes useful to introduce the integral kernel of a density operator 
7 G £?n, which is called a (one-particle) density matrix, and is usually also denoted 
by 7. It is by definition the function 7 G L 2 (R 3 x R 3 ) such that 

V^L 2 (M 3 ), (70)( r )=/ 7(r,r')<KrVr'. (13) 

Jr 3 

The expression of the density matrix 7 in terms of natural orbitals and occupation 
numbers thus reads 

7(r,r') = f «,-</),■« 

i=l 

Formally Py(r) = 7(1", r) and this relation makes sense rigorously as soon as the den- 
sity matrix 7 has a trace on the three-dimensional vector subspace {(r,r),r€l 3 } 
of R 3 x R 3 . 

Let us now clarify the link between the description of electronic structures in 
terms of wavefunctions and the one in terms of density operators. 

The density matrix associated with a wavefunction f G A^[L 2 (R 3 ) such that 
|| f ||z,2( K 3JV) = 1 is the function of L 2 (R 3 x R 3 ) defined as 

7p(r,r')=Jv/ V(r,T 2 ,— ,r N )¥(T f ,T2,— ,r N )dT 2 -dr N (14) 

(recall that we are dealing with real-valued wavefunctions), and the corresponding 
density operator by 

V</>GL 2 (R 3 ), ()V4>)(r) = f Mry)Hr')dr'. (15) 

JR 3 

It is easy to see that the density operator y^ is in £$n. Under the additional assump- 
tion that f G Af =l H l (R 3 ), it is even in Besides, the definition $5^ of the density 
associated with f agrees with the definition $Y2\ of the density associated with yp, 
i.e. 

Pf = Pyy, 

and the same holds with the definition of the kinetic energy if f G /\^ =l H l (R 3 ): 

(«p|r|«p> = rfo F ). 



Remarkl. The maps | f G A£ i^ 2 (R 3 ) | |||f || L 2 (m3W) = 1 j 3 f H> y^ G 2> N and 

I f G /\? =l H 2 (R 3 ) j |||f || L 2 (M 3jv) = lj 9 f ^ yp G are not surjective. This 
means that an element of S^n (resp. of is not necessarily the density operator 
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associated with some pure state. However any 7 g @ N (resp. any 7 6 is the (one- 
particle) density operator associated with some mixed state (represented by some 
Af-particle density operator). This property is referred to as the A^-representability 
property of density operators. 

We can now reformulate the electronic structure problem for a system ofN non- 
interacting electrons, in terms of density operators: 

1. The energy of a wavefunction f e /\f =l H l (R 3 ) is a linear form with respect to 
the density operator 7^: 

(«P|f#|S») = E° pmc fa) where E° pmc (7) = Tr (-^7) + j p y V nuc ; 

2. The ground state density matrix, that is the density operator associated with the 
ground state wavefunction f defined by (0, is the orthogonal projector (for the 
L 2 inner product) on the space Span(0i , • • • , 0jv): 

N 

ft* = Lift) 

3. The ground state energy and the ground state density operators are obtained by 
solving the minimization problem 

inf{^nuo(7), 7 e J^(L 2 (IR 3 )), < y< 1, Tr(y) = Tr(-/ly) < ~J . (16) 

The advantages of the density operator formulation, which are not obvious for finite 
systems, will clearly appear in Section[3] where we deal with crystals. 



2.4 The Hartree model and other density operator models of 
electronic structures 

Let us now reintroduce the Coulomb interaction between electrons, taking as a start- 
ing point the non-interacting system introduced in Section 12.21 The models pre- 
sented in this section are density operator models in the sense that the ground state 
energy and density are obtained by minimizing some explicit functional E p nne(y) 
over the set of A^-representable density operators 

All these models share the same mathematical structure. They read: 
mf{£ p nu C (7), 7 e ^(L 2 (K 3 )), < 7< 1, Tr(y) = N, Tr(-Ay) < <»} , (17) 
with 

Epnuc (7) = Tr (-^Ar) + jTj PrVp"* + \ D i.Py^Py)+E{7), 
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where 

*>(/>*)= f f (18) 
Js?Jm? |r-r'| 

is the classical Coulomb interaction and E(y) some correction term. Note that 
D(f,g) is well defined for / and g in L 6 / 5 (M 3 ), see for instance Ej] Section IX.4]. 
Recall also that for each 7 G &> N , p r G L 1 (R 3 ) HL 3 (R 3 ) L 6 / 5 (R 3 ). 

The Hartree model, on which we will focus in this proceeding, corresponds to 

%)=0: 

E^CyHTrf-W) + f^p y V p ^+ l -D{p y ,p y ). 

The reason why we study this model is that it has much nicer mathematical proper- 
ties than other models with E(y) ^ (see below). 

The Kohn-Sham models [24] originate from the Density Functional Theory 
(DFT) lfl3l . In this kind of models, E(y) is an explicit functional of the density 
py, called the exchange-correlation functional: 

EfLffl^f—Ay) + j^V p ^p y + l -D{p v p 7 )+E^{p y ). (19) 

If follows from the Hohenberg-Kohn theorem ETIl (see 11271 for a more mathemat- 
ical presentation of this result) that there exists some functional E xc (p) depend- 
ing only on the density p, such that minimizing ( TPTl i with E p mic = £^ c provides 
the exact ground state energy and density, whatever the nuclear charge distribution 
pnuc N ote however, that the Kohn-Sham ground state density operator obtained 
by minimizing ( fTTb is not the ground state density operator corresponding to the 
ground state wavefunction f °. Unfortunately, the exact exchange-correlation func- 
tional is not known. Many approximate functionals have been proposed, and new 
ones come up on a regular basis. For the sake of illustration, the simplest approxi- 
mate exchange-correlation functional (but clearly not the best one) is the so-called 
Xa functional 

E£ a (p) = -C Xa [ p"\ 
J®? 

where Cxa is a positive constant 

Lastly, the models issued from the Density-Matrix Functional Theory (DMFT) 
involve functionals £(7) depending explicitly on the density operator 7, but not only 
on the density p y . Similar to DFT, there exists an exact (but unknown) functional 
E(y) for which minimizing ( TT7T > gives the exact ground state energy and density, 
whatever the nuclear charge distribution p nuc . However, unlike the exact DFT func- 
tional, the exact DMFT functional also provides the exact ground state density op- 
erator. Several approximate DMFT functionals have been proposed. Note that the 
Hartree-Fock model, which is usually defined as the variational approximation of 
(0 obtained by restricting the minimization set to the set of finite energy Slater 
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determinants, can also be seen as a DMFT functional 

E™c (7) = Tr (- Uy)+ [ p r V p ™ + l -D(p y , p y ) - I f f ffi^ drdr', 

where, as above, y(r,r') denotes the integral kernel of 7. 

The existence of a solution to WIS for a neutral or positively charged system is 



established in [34] for the Hartree model (£ xc = 0), in 11261 for the Hartree-Fock 
model, in ||4) for the Xa and the standard LDA model, and in iT31l for the Miiller 
DMFT functional. 

The key-property allowing for a comprehensive mathematical analysis of the 
bulk limit for the Hartree model is that the ground state density is unique (which 
is not the case for the other models presented in this section). This means that in the 
Hartree framework, all the minimizers to ( fTTI i share the same density. This follows 
from the fact that the ground state Hartree density solves the variational problem 

M j#(p), p>0, VP e//'(M 3 ), f R S= N }> ( 20 ) 

£{p)=F{p) + f pV p ^+ l -D(p,p) 

Js.3 2 



where 



and 



F{p)=inf^Tr(-^AY^j , ye J^(L 2 (R 3 ), 0< y< 1, Tr(y) =N, Tr(-4y) < °°, Py 
As the functional <§{p ) is strictly convex on the convex set 

[p>0, vpe^(R 3 ), J R3 P= N ) 



uniqueness follows. 

The Euler equation for the Hartree model reads 



7° = X>,|fc>(fc|, P°(r) = Pr o(r) = J>|fc(r)| 2 , 

;=1 i=l 



«/ = 1 if £, < £ F , < «j < 1 if £; = £ F , «i = if £,■ > £ F , n > =N > 



H° = --A+V°, 
2 

I -AV°=4n(p nuc -p°) 
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It can be proved that the essential spectrum of the self-adjoint operator H is equal 
to and that, for a neutral or positively charged system, H° has at least N negative 
eigenvalues. The scalar £p, called the Fermi level, can be interpreted as the Lagrange 
multiplier of the constraint Tr(y°) = N. 

Assuming that £n < £jv+i, the ground state density operator y° of the Hartree 
model is unique: It is the orthogonal projector 

f = £\kM\- 

(=1 

In this case, (f2TT > can be rewritten under the more compact form 
'f = l { ^ M (H% p°=p f , 

< H° = ~A+V°, (22) 

for any < £p < £n+i- In this equation, the notation I £f -^(H ) is used for the 
spectral projector of H° corresponding to the spectrum in the interval (— °°, £p]. 

Lastly, we remark that if smeared nuclei are used, then D(p™f : p™f) is well 
defined (and finite). This allows us to reformulate the Hartree ground state problem 
as 

inf{^? u r e (r), ye .y(L 2 (R 3 )), < y< 1, Tr(y) =N, Tt(-Ay) < ~} , (23) 
where 

Ep™(y) = Tr i^-^y) + \d{p™ c - P r ,P nuc - Pr)- 

The main interest of this new formulation of the Hartree problem is that the func- 
tional £ , pn I S ree is the sum of two non-negative contributions: the kinetic energy and 
the Coulomb energy of the total charge distribution p nuc — Py. The presence of the 
unphysical terms corresponding to the self-interaction of nuclei in D(p™f,p™f) is 
not a problem for our purpose. 

The time-dependent version of the Hartree model formally reads 

-±A-(p™(t)- Pr{t) )*\-\-\y(t) , 

where [A,B] = AB — BA denotes the commutator of the operators A and B. We are 
not going to elaborate further on the precise mathematical meaning of this formal 
equation for finite systems, but refer the reader to [5] and references therein (see in 
particular [ 12, Section XVII. B. 5]) for further precision on the mathematical mean- 
ing of the above equation. On the other hand, we will define and study a mild version 
of it in the case of crystals with defects in Section l4~4l 
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3 The Hartree model for crystals 



The Hartree model presented in the previous section describes a finite system of N 
electrons in the electrostatic potential created by a nuclear density of charge p nuc . 
Our goal is to describe an infinite crystalline material obtained in the bulk limit. 
In fact we shall consider two such systems. The first one is the periodic crystal 
obtained when, in the bulk limit, the nuclear density approaches the periodic nuclear 
distribution of the perfect crystal: 

p nUC ->Pp n eT, (24) 

Ppei C being a ^-periodic distribution. The set S% is a periodic lattice of R 3 : 

^ = Zai+Za 2 + Za 3 , (25) 

where (ai , a 2 , a 3 ) is a given triplet of linearly independent vectors of R 3 . The second 
system is the previous crystal in the presence of a local defect: 

~nuc . ,,nuc i /<->£\ 

P ^Pper+ OT > ( 26 ) 

m representing the nuclear charge of the defect. The functional spaces in which p™ c 
and m are chosen are made precise below. 



3.1 Basics of Fourier and Bloch-Floquet theories 

A perfect crystal is characterized by a lattice 3£ of R 3 and a ^-periodic nu- 
clear charge distribution Pp" r c . Not surprisingly, Fourier and Bloch-Floquet theories, 
which allow to conveniently exploit the periodicity of the problem, play essential 
roles in the mathematical description of the electronic structure of crystals. 

Let St 1 be the reciprocal lattice of the lattice ffi defined in (l25T l (also called dual 
lattice): 

= 1a\ + Za* 2 + Zal , where a,- • a* = 2ffSy . 

Denote by F a unit cell of Recall that a unit cell is a semi-open bounded polytope 
of R 3 such that the cells F + R = {(r + R),reF}forRe^ form a tessellation 
of the space R 3 (i.e. (F + R) n (F + R') = if R ^ R' and U RG .«(F + R) = R 3 ). 
A possible choice for F is {x\ai +X2&2 + ^3^3, — 1/2 < Xj < 1 /2}. Another choice 
is the Wigner-Seitz cell of M, which is by definition the semi-open Voronoi cell 
of the origin for the lattice £%. Lastly, we denote by F* the first Brillouin zone, 
that is the Wigner-Seitz cell of the dual lattice. Let us illustrate these concepts on 
the simplest example, the cubic lattice, for which 8% = ai? (for some a > 0). In 
this particular case, £%* — ^jl? , the Wigner-Seitz cell is F = (—a/2, a/2] 3 and 
F* = (— 7r/a, 7r/a] 3 . 
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For each K £ we denote by eK( r ) = |-T|~ 1//2 e' K ' r the Fourier mode with 
wavevector K. According to the theory of Fourier series, each ^"-periodic distribu- 
tion v can be expanded in Fourier series as 

v= E c K(v)eK, (27) 

Kg** 

where ck(v) is the K-th Fourier coefficient of v, the convergence of the series hold- 
ing in the distributional sense. We introduce the usual .^-periodic LP spaces defined 
by 

i&rCO := { v £ Lf 0C (R 3 ) | v ^-periodic} , 
and endow them with the norms 

ll v llz| er (r) : = (J r Vf^j for 1 < p < oo and \\v\\ L ~ a (r) ■= ess-sup |v|. 
In particular, 

\Hi^{r) = ( v > v )Ll {n where ( v ' w )^(r) : = f/ w - 

Any function v £ L 2 er (.T) can be expanded in Fourier modes according to d27l ). the 
Fourier coefficients being given by the simple formula 

c K (v) = ^ p | r V(r)e-' K ^r, 
and the convergence of the series (|27| | also holds in L 2 er (T). Besides, 



v(v,w) ei^ir) xL 2 v jr), (v,w) L i {r) = £ c K (v)c K (w). 

Kg** 

For each set, the ^"-periodic Sobolev space of index s is defined as 



[ Kg** 



and endowed with the inner product 



I (i+iK| 2 r| CK (v)i 2 <. 

kg** 



(v.wj^r) := £ (1 + |K| 2 ) s Ck (v)ck(w). 

Kg** 



The Bloch-Floquet theory was introduced by Floquet for periodic differential 
equations and generalized by Bloch to periodic partial differential equations. We 
just recall the basic results of this theory used in this proceeding and refer the reader 
to |3 1 1 for further precisions. 

Any function / £ L 2 (R 3 ) can be decomposed by the Bloch-Floquet transform as 



16 Eric Cances, Mathieu Lewin and Gabriel Stoltz 

/(r)=^/ q (r)e'' qr dq, 

where j T , is a notation for \T* |~ 1 J r , and where the functions / q are defined by 

Mr) = E /(r + R)e-"»-("-+ R ) = £ /(q + K)e^. (28) 

re,* I 1 I Kg^* 

For almost all q G K 3 , / q G L 2 er (r). Besides, / q+K (r) = / q (r)e _iK-r for all Kef* 
and almost all q G K 3 . Lastly, 

For R G K 3 , we denote by Tr the translation operator defined by 
VveL 2 (R 3 ), (T R v)(r)=v(r-R). 



The main interest of the Bloch-Floquet transform ( 1281 ) is that it provides a "block di- 
agonalization" of any ^-periodic operator, that is of any operator on L 2 (IR 3 ) which 
commutes with Tr for all Rgf. Consider first a bounded ^-periodic operator A 
on L 2 (R 3 ). Then there exists a family (A q ) qe r* of bounded operators on L 2 el .(.T) 
such that 

VveL 2 (R 3 ), (Av) q =A q Vq for almost all q G T* . (29) 

If, in addition, A is self-adjoint on L 2 (R 3 ), then A q is self-adjoint on L 2 er (F) for 
almost all q G r* and 



<T(A) = U ^(A,). 

qer* 

In particular, the translation operators (Tr)r £ ^, which obviously commute with 
each other, are homotheties in the Bloch-Floquet representation 

VRGf, (T R ) q = e^ R l z|er(r) . 

As (^k)kg^* form an orthonormal basis of L 2 er (F), it follows from (129b that any 
bounded ^-periodic operator on L 2 (R 3 ) is completely characterized by the Bloch- 
Floquet matrices (([A K ,k'(q)])(k .k')g,**x ,??*)qer* defined for almost all q G T* by 

^K,K'(q) := (<?K,Aqe K ')L2 cr (r)- 

In particular, it holds 

VveL 2 (K 3 ), (A7)(q + K)= £ A K , K ,(q)v(q + K'), 

for all (K,K') ef'xf and almost all qef*. 
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For unbounded operators, the situation is a little bit more intricate. Let us limit 
ourselves to the case of^-periodic Schrodinger operators of the form 

H = ~A+V pel 

with Vper € ip er (r). By the Kato-Rellich theorem and ED Theorem XIII.96], the 
operator H is self-adjoint on L 2 (R 3 ), with domain // 2 (M 3 ). It can also be decom- 
posed as follows: 

VvG// 2 (K 3 ), v q £ Hpe r (r) and (//v) q =H q v q for almost all qeT, 

where i/ q is the self-adjoint operator on L 2 er (T) with domain i/ 2 el .(r), defined by 

1 Iql 2 
H q = --A-iq-V + ^-+V peI . 

It is easily seen that for each q e r*, H q is bounded below and has a compact 
resolvent. Consequently, there exists a sequence (e„. q )«>i of real numbers going to 
+°o, and an orthonormal basis (u„ jq ) n >i of L 2 er (T) such that 

+<*> 

= ^ £«, q | M «, q ) ( M n.q| • 
n=\ 

As the mapping q n- // q is polynomial on R 3 , it is possible to number the eigenval- 
ues £„ q in such a way that (e„.o) n >i is non-decreasing and that for each n > 1, the 
mapping q H> e„. q is analytic in each direction. Then (see Fig.[TJ) 

a(H) = jj a(H q ) = (J , 
q er* «>i 

with 

£ « = m in e«.q, ^ = max£„ iq . 
q er* qer* 

The interval is called the n th band of the spectrum of H 

prove that the spectrum of H is purely absolutely continuous ll35l 
has no eigenvalues. 



(30) 

. It is possible to 
. In particular, H 



3.2 Perfect crystals 

The purpose of this section is to formally construct, then justify with mathematical 
arguments, a Hartree model for the electronic structure of perfect crystals. 

As announced, we begin with a formal argument and consider a sequence of 
finite nuclear distribution (p" uc )neN converging to the periodic distribution p^f of 
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Fig. 1 The spectrum of a periodic Schrodinger operator is a union of bands, as a consequence of 
the Bloch-Floquet decomposition. 



the perfect crystal when n goes to infinity. For instance, we can take 



Pn 




(we assume that the function describing the nuclear charge in the unit cell of the 
perfect crystal is supported in some compact set included in the interior of F). We 
solve the Hartree problem for each p" uc with the constraint that the system remains 
neutral for each n. Assuming that when n goes to infinity, 



the Hartree ground state density converges to some ^-periodic density Pp er G 



• the Coulomb potential generated by the total charge converges to some 
periodic potential Vp er ; 

• the Hartree ground state density operator converges to some operator }£ er ; 

• the Fermi level converges to some £p 6 R, 

we obtain by formally passing to the limit in d22L the self-consistent equations 

T/per = 1 (-«,, £ 0] (#per) , Pper = P$ er , 



-AV° ei = 4n(p- 



1 



V 



per - 

nuc _ 
per 



(31) 



Pp°e r )- 



Let us comment on this system of equations. First, we notice that for the periodic 
Coulomb equation — AV® a = 47f(p 
be neutral: 



^nuc 
J per 



Pper) t0 nave a solution, each unit cell must 



o nuc - Z 

Pper — 



(32) 
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where Z is the number of electrons, and also the number of protons, per unit cell. 
Second, as Vp er is ^-periodic (and belongs to Lp el .(F) even for point-like nuclei), 
we can apply the result of the previous section and write down the Bloch-Floquet 
decomposition of H^ er : 

«r)q = - ^-jq-V + l^+y p ° er = £e n , q Kq)K q |. (33) 

n=l 

The operator y^ el = £ oj (H^ el ) then is a bounded self-adjoint operator which 
commutes with the translations (Tr)r & ^>, and its Bloch-Floquet decomposition 
reads 

(^er)q= L 1 e„, q <£° Kl) Kll ' 
n=\ 

Actually, the set {q g F* \ 3n > 1 s.t. £„ q = £p} is of measure zero (the spectrum 
of //p er is purely continuous). It follows that y^ er is always an orthogonal projector, 
even if Ep 1 belongs to the spectrum of //p er . 

Using the Bloch decomposition of ^ er , we can write the density Pp er as 

Pper(r) = / E l e n . q <s»\ u >M\ 2 dq- 
Jr* n =\ 

Integrating on F, and using d32l > and the orthonormality of the functions (w M .q)«>i 
in Lp er (F), we obtain 

Z =TF^EK < l er * I £ «-^ £ f}|- (34) 

I 1 I n=\ 

This equation determines the value of the Fermi level £p uniquely. It is easy to see 
that if the periodic Coulomb potential is shifted by a uniform constant C, and if gp 1 
is replaced with £p + C, then and Pp er remain unchanged. 

The formal bulk limit argument presented above has been rigorously founded by 
Catto, Le Bris and Lions in lITTl . for p™ c = Lrgz 3 X(' ~ R) (smeared nuclei of unit 
charge disposed on the cubic lattice Z 3 ). It is also possible to justify the periodic 
Hartree model by passing to the limit on the supercell model with artificial periodic 
boundary conditions (see [7|). The latter approach is less physical, but technically 
much easier, and its extension to arbitrary crystalline structures (including point-like 
nuclei) is straightforward. It results from these mathematical works that the Hartree 
model for perfect crystals is well-defined. More precisely: 

1 . The Hartree ground state density operator )^ er and density Pp er of a crystal with 
periodic nuclear density p™ r c (composed of point-like or smeared nuclei) are 
uniquely defined; 

2. The ground state density Pp el . satisfies the neutrality charge constraint d32l ; 
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3. The periodic Coulomb potential V^ r and the Fermi level £p are uniquely defined 
up to an additive constant (and Vp er — £p is uniquely defined); 

4. The ground state density operator y^ eT is an infinite rank orthogonal projector 
satisfying the self-consistent equation (l3TT i; 

5. y® er . can be obtained by minimizing some periodic model set on the unit cell F 
(see HB for details). 

In the remainder of the paper we assume that the system is an insulator (or a 
semi-conductor) in the sense that the N th band is strictly below the (A^+ l) st band: 

E + < E~ 

where E^ are defined in (f30b - In this case, one can choose for £p any number in the 
range (E^ ,E^ +l ). The electronic state of the perfect crystal is the same whatever 
the value of £p in the gap (E^ ,E !V+l ). On the other hand, as will be seen in the next 
section, fixing the value of £p may change the electronic state of the crystal in the 
presence of a local defect. 

In this paper however, we are only interested in the dielectric response of the 
crystal, which corresponds to the limit of small defects (in a sense that will be made 
precise later), and in this limit, the value of £p does not play any role as long as it 
remains inside the gap (E^ ,E^ +l ). For simplicity, we consider in the following 

f + -+- y~ 

%- j ' 

Lastly, we denote by 

g = E N+l -E+>Q (35) 

the band gap. 



3.3 Crystals with local defects 

We now describe the results of Q dealing with the modelling of local defects in 
crystals in the framework of the Hartree model. The main idea is to seek the ground 
state density operator of a crystal with a local defect characterized by the nuclear 
charge distribution ( T26T ) under the form 




In this formalism, the defect is seen as a quasi-molecule with nuclear charge dis- 
tribution m and electronic ground state density operator Q £ o (and ground state 
electronic density On „ ), embedded in the perfect crystal. Here, the charge of the 

defect is controlled by the Fermi level (the chemical potential). The dual approach, 
in which the charge of the defect is imposed, is also dealt with in J7J. It should be 
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noticed that neither m nor On n are a priori non-negative. For instance, the nuclear 

ffl,Sp 

distribution of a defect corresponding to the replacement of a nuclear of charge z 
located at point R € R 3 with a nucleus of charge z' is m = (z' — z)5r and can there- 
fore be positively or negatively charged depending on the value of { — z. Regarding 
the electronic state, the constraints (y m£ o)* = y m6 o,0<y me o<l and p r . > 0, 

' F ' F ' F fl!i£p 

respectively read (Q m e oJ = Q m & -$, r < Q m £ a < 1 - ^ and p Q > -p° er 



The next step is to exhibit a variational model allowing to compute Q o from 

m, £p and the ground state of the perfect crystal. 

First, we perform the following formal calculation of the difference between the 
Hartree free energy of some trial density operator y = ^ er + Q subjected to the 
nuclear potential generated by p™ r c + m, and the Hartree free energy of the perfect 
crystal: 

(EpgTA + Q) - + Q)) (^g-O&r) - 



fonnal 



Tr (-~AQ) + p e V p er - f p Q V m + l -D{p Ql p Q ) - e°Tr(g) 

mV° er +~Z)(m,m). (36) 

The last two terms are constants that we can discard. Of course, the left-hand side of 
(f36l > does not have any mathematical sense since it is the difference of two energies 
both equal to plus infinity. On the other hand, we are going to see that it is possible 
to give a mathematical meaning to the sum of the first five terms of the right-hand 
side when Q belongs to some functional space =2 defined below, and to characterize 
the ground state density operator Q £ o of the quasi-molecule, by minimizing the 
so-defined energy functional on a closed convex subset ,Jf of =2. 

For this purpose, we first need to extend the definition ( fT8b of the Coulomb in- 
teraction to the Coulomb space ^ defined as 



feLl M 3 ),D(fJ):=4nJJ-^dk} 



where y'(R 3 ) is the space of tempered distributions on R 3 . Endowed with its nat- 
ural inner product 



(f,gW :=D(f,g) :=47t [ ^§^dk, (37) 

^ is a Hilbert space. It can be proved that L 6 / 5 (R 3 ) <-t and that for any (f,g) € 
L 6 / 5 (R 3 ) xL 6 / 5 (R 3 ), it holds 



Jr3 W J®?Jm? |r-r| 
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Hence, the definition (f3Tb of D(-, •) on is consistent with the usual definition ( fTST l 
of the Coulomb interaction when the latter makes sense. The Coulomb space ^€ 
therefore is the set of charge distributions of finite Coulomb energy. 
Second, we introduce, for an operator A on L 2 (M 3 ), the notation 

A-:=^ er A^ a , A-+:=^ er A(l-^ er ), 
A + - := (1 - tfjA^, A++ := (1 - )pA(l - )p, 

and note that the constraints Q = Q* and — }^ er < Q < 1 — )^ el . are equivalent to 

Q* = Q, Q 2 <Q ++ ~Q-- (38) 

From the second inequality we deduce that it then holds Q < and Q ++ > 0. 
Using the fact that Tr(V p ° el .g) = J r3 PqV^ si , we formally obtain 

Tr (-\aq\ + j p e y» r - 4 } Tr(e) = Tr((//» er - eg)fi) 

= Tr((// p ° er - e F °) ++ e ++ ) +Tr((// p ° er - $-(2-). 

We now remark that, by definition of $, r , {H° ei ~ e°)++ > and (H^-eg)— <0, 
so that the right-hand term of the above expression can be rewritten as 

Tr(|< r - eF | 1/2 (e ++ -e-)l< r -e F °| 1/2 )- (39) 

The above expression is well defined in M + U {+°°} for all Q satisfying the con- 
straints ( |38l ). It takes a finite value if Q is chosen in the vector space 

^ = {Qe& 2 \Q* = Q, Q—eB u Q ++ eei, (40) 
|V|ee6 2 , |v|fi— |V|eSi, |V|e++|V|e6i}, 

where ©i and ©2 respectively denote the spaces of trace-class and Hilbert-Schmidt 
operators on L 2 (M 3 ) (see Appendix for details). Endowed with its natural norm, or 
with any equivalent norm such as 

ll<2iu=ll(i+|v|)ei| e2 +ii(i+|v|)e ++ (i+|v|)i| ei +ii(i+|v|)(2-(i+|v|)i| ei , 

£2, is a Banach space. 

Before proceeding further, let us comment on the definition of £2. As the trial 
density operators Q must satisfy the constraints (1381 . it is natural to impose Q* = Q. 
Since |//p er — £p| 1//2 (l + |V|) _1 is a bounded operator with bounded inverse (see Q), 
the four conditions Q~ e 61, Q++ G 61, |V|Q~ |V| e ©j and |V|g++|V| e ©[ 
are necessary and sufficient conditions for the expression ( f39b with Q satisfying d38l 
being finite. The other constraints imposed to the elements of i? (that is, Q £ 62 
and I V\Q g ©2) follow from the fact that for any Q satisfying (l38l 
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(Q~~ e ©i, Q ++ e 61) => (e 2 e©i) 
(]v|Q— ]V| G ©1, |V|G ++ |V|G6i) (|V|G 2 |V| G 61) . 

In order to simplify the notation, we set for Q E J2, 

Tr (e):=Tr(e++ + e-), 

Tr (« r - e F °)e) := Tr(|H p ° er - £ F °| 1/2 (e ++ - Q~)^ %T /2 )- 

An important result is that the linear application Q 1— > Pq originally defined on the 
dense subset J2n&i of £i can be extended in a unique way to a continuous linear 
application 

Q >-> Pq- 

Note that the density associated with a generic element of i? is not necessarily an 
integrable function. On the other hand, its Coulomb energy is always finite. 

Let m be such that V m = (m* \ ■ € Here and in the sequel 
:= {y g L 6 {R 3 ) I VV G (L 2 (R 3 )) 3 } 
denotes the dual space of If, endowed with the inner product 

<yi,v 2 )* := f vy 1 -vy 2 = -i- f \k\ 2 f^k)v 2 (k)dk. 

4n Jr3 47T J R 3 

It follows from the above arguments that the energy functional 

E m ^(Q) = Tr (« r - e F °)e) - / p Q V m + \d{ PQiPq ) 

is well defined on and that a good candidate for a variational model allowing to 
compute the ground state density operator Q m e o is 

M{F"^(<2),<2GJf} (41) 

where 

jr={Q€4\-'£ t <Q<l-'& t }. (42) 

Note that is a closed convex subset of J2. 

The above formal construction of the model (RTt is justified in [7] by means of 
rigorous bulk limit arguments. To summarize the situation, the Hartree ground state 
density operator of the crystal with nuclear charge density Pp" r c + m (the charge of 
the defect being controlled by the Fermi level) is given by 
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Y=1@* + Qm4 

where Q £ o is obtained by solving (|4TT >. 

The existence of a Hartree ground state density operator for a crystal with a 
local defect, as well as the uniqueness of the corresponding density and some other 
important properties, are granted by the following theorem which gathers several 
results from Q and |9l . 

Theorem 1. Let m such that (m* | • G L 2 (M?) + Then, 

1. ( 14 li has at least one minimizer Q m o, and all the minimizers of (141\) share the 
same density p „o," 

2. Q m gO is solution to the self-consistent equation 

Qm4 = (<r + (P m , £ o • I" 1 ) - «r) +5, (43) 

where 8 is a finite-rank self-adjoint operator on L 2 (R 3 ) such thatO < 8 < 1 and 
Ran{8) C Ker (i/ p ° er + (p m £ o - m) * | • |" 1 - eg) . 

The interpretation of the Euler equation (l43l , which also reads 
with 

K4 = <r + (P m , e o -m)*| • I" 1 , < 5 < 1, Ran(5) C Ker(^ - fig), 
is the following. The mean-field Hamiltonian H° is uniquely defined, since all 

m,£p 

the minimizers of < f4TT> share the same density p m £ o . Besides, the operator (p m £ o — 
m) ★ I • |~' being a relatively compact perturbation of H® el , it results from the Weyl 
theorem (see ll3~Tl Section XIII. 4]) that the Hamiltonians and H° have the 

same essential spectra. On the other hand, while H® er has no eigenvalues, H° £0 may 

have a countable number of isolated eigenvalues of finite multiplicities in the gaps as 
well as below the bottom of the essential spectrum. The only possible accumulation 
points of these eigenvalues are the edges of the bands. 

If eg ^ go)> men ^ = and the ground state density operator of the crystal 

in the presence of the defect is the orthogonal projector }^ er + Q m £ o: All the energy 
levels lower that the Fermi level are fully occupied while the other ones are empty 
(see Fig. O. In this case, Q £ o is both a Hilbert-Schmidt operator and the difference 
of two projectors. It therefore follows from [18, Lemma 2] that 



Tro(G ffl , £ o)eN. 



(44) 
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Assuming that m e L l (M?) and J r3 m € N, the integer 




can be interpreted as the bare charge of the defect (in contrast with the screened or 
renormalized charge to be defined later). 

If £° e ), then the energy levels with energy ep* may be fully or partially 

occupied, and it may a priori happen that d4TT > has several minimizers, differing from 
one another by a finite rank self-adjoint operator with range in Ker(i/° — eS). 

fc F 

^ 1 h-> -«, E o) 

I | 1 1-rH 

Fermi sea electrons 

Fig. 2 General form of the spectrum of the self-consistent operator ^ , in the presence of a 
defect and for a fixed chemical potential Bp. 



4 Dielectric response of a crystal 

In this section, we study the response of the electronic ground state of a crystal to a 
small, effective potential. In Section |4~T1 we consider a time-independent perturba- 
tion V <EL 2 (R 3 ) + % J ' , with 1 1 V | L 2 +% f/ < a (for some a > small enough). It can be 
proved (see [9 , Lemma 5]) that there exists j3 > such that 

(ihM-r'w^Kp) => (ii(p B ,^-»)*i-r 1 iii*+»'< a )- ( 45 ) 

The results of Section 14.11 therefore directly apply to the case of a crystal with a 
local defect with nuclear charge distribution m, provided the defect is small enough 
(in the sense that \\m~k \ ■ | _1 ||i2 + ^/ < /3). 

In Section l4~4l we consider a time-dependent perturbation 

v(f,r) = (p(f,-)*|-r 1 )(r) with pei(l,L 2 (R 3 )nf). (46) 



4.1 Series expansion of the time-independent response 

For V € L 2 (R 3 )+ c io' , the spectrum of H® m + V depends continuously of V. In par- 
ticular (see [9 , Lemma 2]), there exists some a > 0, such that if £ is a smooth curve 
in the complex plane enclosing the whole spectrum of H® er below gp, crossing the 
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real line at ep 1 and at some c < infer (i/p er ) and such that 

d(a{H^),A) = | where A = [z e C | d(z, t) < | } , 



d denoting the Euclidean distance in the complex plane and g the band gap (f35T > 
(see Fig. |3J, then (j{H® er +V) n (— °°,£p] is contained in the interior of <£ for all 



V SL 2 (R 3 )+^" such that H^IU^/ < a. 




As a consequence, we obtain that for all V € L (R ) + < ^" such that ||V||£2 + <g>/ < a 



per/ 



2/7T 



(47) 



where we have used the fact that £p ^ cr(i/p er + V) to establish the first equality, and 
the Cauchy formula to derive the second one. 
Expanding d47b in powers of V, we obtain 



Qv = Qn.V + QN+l,V, 
n=\ 



(48) 



where we have gathered the terms involving powers of V larger than N in a remain- 
der Qn+i.v- The linear contribution is given by 



1 



Qlv = ^f (z-H^) Wd-HL) 1 dz. 



per/ ' V per/ 

The higher order contributions and the remainder are respectively given by 



and 



- 2^ V) ' 



v z-h; 



dz 



N+l 



per/ 



(49) 
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Proposition 1. The terms of the perturbation expansion ( 1451 ) enjoy the following 
properties. 

1. The k-linear application 

(Vl , • • • ,V n ) H> ^ £ (z - // p ° er ) ~ ' Vi (z - ff p ° er ) ~ ' ■ ■ ■ V„ (z - < r ) " 1 dz 

is well-defined and continuous from (L 2 (R 3 ) + e £") n to J2 for all n > 1, and from 
(L 2 (R 3 )+f')"to6i /or a« n>6. In particular, for all V € L 2 (R 3 ) + 'if, Q„,y G 
i?/or alln>\ and Q„y £ ©i /or all n>6. Besides, for all V € L 2 (R 3 ) + <£' , 
Tro(2n,v) = for all n > 1 and Tx{Q„y) = for all n>6. 

2. IfV € L 1 (R 3 ), g n ,y is in &\for each n>\and Tx(Q n y) = 0. 

3. For each V €E L 2 (K 3 ) such that \\V\\ L 2 . < G!, f/ze operator Qn+i.v is i n & 
for allN > vw'f/z Tro(j2iv+i.v) — 0, and in &\ for all N > 5, with Tr(gAr+i,y) = 

Ttt)0av + i,v)=o. 

We are now in position to define some operators which play an important role in 
the sequel: 

• the Coulomb operator v c , which defines a bijective isometry between and c €'\ 

Vc (p) : = p*|-r'; 

• the independent particle polarization operator Xo defined by 

Xo(V) :=pQ ltV , 

which provides the first order response of the electronic density of the crystal to 
a time-independent modification of the effective potential. The operator Xo is a 
continuous linear application from L*(R 3 ) to L^R 3 ) and from L 2 (R 3 ) + e £" to 
L 2 (R 3 )n^; 

• the linear operator Jzf defined by 

■&■■= -Xovc 

which is a bounded nonnegative self-adjoint operator on c <o . As a consequence, 
( 1 + _Sf ) ~ 1 is a well-defined bounded self-adjoint operator on ^ ; 

• the dielectric operator £ = v c (l +^f)v ( T 1 , and its inverse, the dielectric permit- 
tivity operator 

e- l =v c {\+X)- l v-\ 
both being continuous linear operators on & . Note that the hermitian dielectric 

~ —1/2 1/2 

operator, defined as e = v c ' ev c is a self-adjoint, invertible, bounded operator 
on L 2 (R 3 ) and is for this reason conveniently used in mathematical proofs. 

We now focus our attention on the total Coulomb potential 

V m = (m - p m£< > )*|.|- 1 = Vc (m- p m£ o ) , 
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generated by some charge distribution m such that 1 1 m * | • | ~ 1 1 1 L 2 < j3 , and on the 
response p £ o of the Fermi sea. In view of d45t . we can apply the above results and 
deduce from fel l that 

Pm4 = PQ-v m =PQu-v,„ +PQ 2 ^ Vm =-X0V m + PQ 2 _ Vm 

= ^ m -P,^ +P Q,-v,„- W 
The above relation, which also reads 

(m- Pm£ o) = (l+J?)- 1 m-(l+J?)- l ( P Q ) (51) 

or 

V m = v c {\+X)- l m-v c {\+X)-\p d ), (52) 

is fundamental since it allows to split the quantities of interest (the total charge 
(m — p m £ o ) or the total Coulomb potential V m generated by the defect) into two 
components: 

• a linear contribution in m, very singular, and responsible for charge renormaliza- 
tion at the microscopic level, and for the dielectric properties of the crystal at the 
macroscopic level; 

• a nonlinear contribution which, in the regime under study (||m*| • | || ji^gi < /3 ), 
is regular at the microscopic level and vanishes in the macroscopic limit. 



4.2 Properties of Q m P o and p m „o for small amplitude defects 



The relation ( l50t , combined with the properties of the operator «5f stated in Propo- 
sition |2] below, allows to derive some interesting properties of Q £ o and p £ o and 
to propose a definition of the renormalized charge of the defect. 

Proposition 2. Let p e L'(R 3 ). Then, Jz? (p) G L 2 (R 3 ) ntf, jSf(p) is continuous on 
R 3 \ ffl*, and for all a € S 2 (the unit sphere ofR 3 ), 

lim J?(p)(j]0) = (a T La)p{0) (53) 
where L € R 3x3 is the non-negative symmetric matrix defined by 

2 



% N 



VkeR 3 , k T m=^-Z £ 



((k- V r )w„ iq ,M M / iq )) L 2 (r) 



^"1 n=ln'=N+\^ r * ( £ n',q~ £ «,< 



3 



-^q, (54) 



where the £„. q 's and the u nSi 's are the eigenvalues and eigenvectors arising in the 
spectral decomposition d33l l of(H^ er ) q . Additionally, 
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L = ~Tr(L) > 0. (55) 

Notice that the convergence of the series ( 154-b is granted by the fact that e n < „ — 
e„, q >E",—L^>g for all n < N < n' and all qef* (where g > is the band gap), 
and the existence of C € R + such that ||« n ,q||//2 ^ < C for all 1 < n < N and all 

qgT*. Actually, the convergence of the series is rather fast since Z7 ~ Cn' 2 ^ 

n'— 7-°° 

(this estimate is obtained by comparing the eigenvalues of H® el with those of the 
Laplace operator on Lp el .(F)). 

We do not reproduce here the quite technical proof of Proposition^ Let us how- 
ever emphasize the essential role played by the long range character of the Coulomb 
potential. If | • | _1 is replaced by a potential v r G L'(R 3 ), then for all p <E /^(R 3 ), 
p * v r e L 1 (R 3 ), hence if (p) € L 1 (R 3 ) and L = 0. More precisely, the Bloch-Floquet 
decomposition of the Coulomb kernel reads 

47T / 1 ^ e' K r \ 



and only the singular component Tfroi! which originates from the long-range of 
the Coulomb potential, gives a nonzero contribution to L. 

We can deduce from (150b and Propositionl2lthat, in general, the minimizer Q e o 
to (|4H is not trace-class and that the density p £ o is not an integrable function if 
the host crystal is anisotropic. Let us detail this point. 

Consider some m e L 1 (R 3 ) nL 2 (M 3 ) such that / R3 m^0and ||m*| • | _1 ||z,2 + ^/ < 
/3. In view of (PES! and Proposition [T] it holds 

Tr o(G m . £ o) - TroCfii,-^ + Qi,- Vm ) = 0. (56) 

Assume that p £ o is in L 1 (R 3 ). Then a technical lemma (see J9] Lemma 4]) shows 
that the Fourier transform of the density p^ ^ , corresponding to the nonlinear 
response terms, is continuous and vanishes at zero. This means that, although it 
is not known whether p= is in L ! (R 3 ), this density of charge behaves in the 
Fourier space as if it was integrable with an integral equal to zero. It follows from 
( f50b and PropositionQ]that for each a e S 2 , 

P m) e°(0)= lim ^(pIX-m)(Tj(T)-((T r L(T)(p o(0)-m(0)). (57) 

As by assumption m(0) ^ (since L 3 m ^ 0), we reach a contradiction unless the 
matrix L is proportional to the identity matrix. Defining here an isotropic crystal as 
a crystal for which L ^ LqI, this proves that, in general, p me o is not an integrable 
function for anisotropic crystals (and this a fortiori implies that Q o is not trace- 
class). 
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Let us now consider an isotropic crystal. If Q m £ o were trace-class, then p m £ o 
would be in L 1 (R 3 ), and we would deduce from d56l > that 

(2^) 3/2 P m , £ «(0) = P m4 = Tr(Q m £ o) = Tr (e m , £ o) = 0. 

Again, except in the very special case when L = 1, this contradicts d57l i since 
m by assumption. Thus, in general, Q g o is not trace-class, even for isotropic 
crystals. We do not know whether the electronic density p o generated by some 

1 F 

m e L l (R i )nL 2 (R i ) (this assumption implies m e L 6/,5 (IR 3 ) ^ ^) in an isotropic 
crystal is integrable or not. If it is, it follows from ( |57T ) that, still under the assump- 
tion that \\m* | • | _1 || L 2 + ^/ < j8, 

f f Is? m 

This quantity can be interpreted as the renormalized charge of the defect, which dif- 
fers from the bare charge J K 3 m — Tro((2 m £ o) = J R 3 m by a screening factor j-^- 
This is formally similar to the charge renormalization phenomenon observed in 
QED (see [ 17] for a mathematical analysis). 



4.3 Dielectric operator and macroscopic dielectric permittivity 

In this section, we focus again on the total potential 

V m = (m-p m4 )*\-\- 1 (58) 

generated by the total charge of the defect, but we study it in a certain macroscopic 
limit. 

For this purpose, we fix some m e L l (M 3 ) n L 2 (R 3 ) and introduce for all 77 > 
the rescaled density 

m n (r) := r] 3 m(rir). 

We then denote by V,H the total potential generated by m-q and the corresponding 
electronic polarization, i.e. 

V^:=(mr,- Pmn4 )*\-\-\ (59) 
and define the rescaled potential 

W%{r):=r]- l V? n (r}- l v). (60) 
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The scaling parameters have been chosen in a way such that in the absence of di- 
electric response (i.e. for _£f = and pg „ — 0), it holds W„? = v c (m) — m* \ ■ | _1 

for all 77 > 0. To obtain a macroscopic limit, we let 77 go to zero. 

As ||(ffJr/ *| • l -1 )!!^' = \\ m rj\\ c e — ? 7 1 ^ 2 |l m l!'^' we can a PPly tne results of the 
previous sections as soon as rj is small enough. Introducing the family of scaling 
operators (t/jj)tj>o defined by (i/^ /)(r) = ?7 3/ ' 2 /( , 7 r ) (each U-q is a bijective isome- 
try of L 2 (R 3 )), the equation linking the density of charge m to the rescaled potential 
Wm reads 

W% = v\ l2 U*£- l U nV l J 2 m + wl, (61) 

where the nonlinear contribution w m is such that there exists C £ R+ such that for 
77 small enough, ||wm||<g>' < Crj. The macroscopic limit of Wm therefore is governed 
by the linear response term, and is obtained as the limit when 77 goes to zero of the 
family (f/^e~ i/tj)rj>o of bounded self-adjoint operators on L 2 (R 3 ). 

If £ 1 was translation invariant, that is, if it was commuting with all the trans- 
lations Tr for R £ R 3 , it would be a multiplication operator in the Fourier space 

(i.e. such that for all / £ L 2 (R 3 ), (£- 1 /)(k) = ^-'(^/(k) for some function 

R 3 3kh>r'(k) G C). Using the fact that the operator v l J 2 is the multiplication 
operator by (4^)'/ 2 /|k| in the Fourier space, we would obtain in the limit 



lim _ ' W m (k) =47rm(k). 

As the operator £ actually commutes only with the translations of the lattice 
the above argument cannot be applied. On the other hand, it can be proved, using 
Bloch-Floquet decomposition, that Wm has a limit W m when rj goes to zero, and that 
this limits satisfies 



A^AiF^J w " ,(k)=4 ^ (k) ' (62) 

where [e _1 ]oo(q) is the entry of the Bloch matrix of the ^"-periodic operator £ _1 
corresponding to K = K' = 0. Besides, 

lim f r~J k|2 /„.^ - ^Mk, (63) 

TJ- " 



0+ V[e _1 ]oo(J7k) 



where £m is a 3 x 3 symmetric, positive definite matrix. Transforming back ( 1621 in 
the physical space, we obtain the macroscopic Poisson equation (O. Let us formalize 
this central result in a theorem. 



Theorem 2. There exists a 3 x 3 symmetric matrix £m > 1 such that for all m £ 
l} (R 3 )nL 2 (R 3 ), the rescaled potential W,H defined by d6Ql ) converges to W m weakly 
in 'to' when 7] goes to zero, where W m is the unique solution in 'i§' to the elliptic 
equation 
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-div(e M VW,„) =4nm. 

The matrix £m is proportional to the identity matrix if the host crystal has the sym- 
metry of the cube. 

From a physical viewpoint, the matrix £m is the electronic contribution to the 
macroscopic dielectric tensor of the host crystal. Note that the other contribution, 
originating from the displacements of the nuclei [29], is not taken into account in 
this study. 

The matrix £m can be computed from the Bloch-Floquet decomposition of H® er 

—1/2 1/2 

as follows. The operator £ = v c £v c ' being ^-periodic, it can be represented by 
the Bloch matrices ( [£kk' (q)]k,k'g^* )qer* ■ It is proven in [9 1 that each entry of the 
Bloch matrix %,k' C 1 ?* 7 ) h as a limit when 77 goes to + for all fixed a 6 S 2 . Indeed, 

lim Eoo(r)o) = 1 + a T Lo 

77^0+ 

where L is the 3x3 non-negative symmetric matrix defined in d54l ). When K, K' ^ 0, 
£ K,K'( I 7 C7 ) nas a limit at 77 = 0, which is independent of a and which we simply 
denote by £k.k' (0)- When K = but K' 7^ 0, the limit is a linear function of a: for 
allK'e^*\{0}, 

lim eoK'(r}a) = Pk>-o, 

T)->0+ 

for some J3 K ' € C 3 . Both £ KK '(0) (K,K' ^ 0) and )3k can be computed from the 
eigenvalues £„ q and eigenvectors u„ y q of the Bloch-Floquet decomposition of H® er 
by formulae similar to dS4b . As already mentioned, the electronic contribution to the 
macroscopic dielectric permittivity is the 3x3 symmetric tensor defined as [|6) 

VkeM 3 , k r £ M k= lim |k|2 . (64) 

By the Schur complement formula, it holds 

K.KV0 



-^ rTr — — = £ 00 (Tjk)- £ £ , K (i7k)[C(T7k) 1 ] KlB »%/, (i7k) 



where C(T?k) 1 is the inverse of the matrix C(rjk) = [£KK'( r ?k)] K K / & ^*y{ }. This 
leads to 

Ikl 2 - ,2 , , r rI , ^r^-n 



lim i 1 ^|k| 2 + k r Lk- £ (/3K-k)[C(0)- 1 ] KiK ,(/3 K -k) 

where C(0)~ l is the inverse of the matrix C(0) = [£kk'(0)]k,k'g,«*\{0}- Therefore, 
£ M = 1+L- £ /3K[C(0)- 1 ] K .K'j3K'- (65) 

K,K'G,«*\{0} 
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As already noticed in (6), it holds 



1 < e M < 1 +L. 



Formula d65l ) has been used in numerical simulations for estimating the macro- 
scopic dielectric permittivity of real insulators and semiconductors J6] [22] [23] [14J 
[16). Direct methods for evaluating £m, bypassing the inversion of the matrix C(0), 
have also been proposed [32 25]]. 



4.4 Time-dependent response 

We study in this section the variation of the electronic state of the crystal when the 
mean-field Hamiltonian H® CT of the perfect crystal is perturbed by a time-dependent 
effective potential v(f,r) of the form d46b . The mathematical proofs of the results 
announced in this section will be given in ifTOl . 



Under the assumption that Pp" r c g Lp er (F) (smeared nuclei), the mean-field potential 
y per is ^-periodic and in C°(R 3 ) nL°°(R 3 ). Besides, there exists a constant C > 
such that ||p * | ■ < C||p|| L 2 n ^ for all p € L 2 (R 3 ) n c i, so that the time- 

dependent perturbation v belongs to L* (R,L°°(M 3 )). 

Let us now define the propagator (^(^^(j^eKxR associated with the time- 
dependent Hamiltonian H v {t) following lf30l Section X.12]. To this end, consider 
first the propagator Uo(t) = e^ ,fH P er associated with the time-independent Hamilto- 
nian //£,,., and the perturbation in the so-called interaction picture: 



Standard techniques (see for instance l28l Section 5.1]) allow to show the existence 
and uniqueness of the family of unitary propagators (i/i n t(f,A'))(.M)eKxH associated 
with the bounded operators (v; nt (f )) fG R, with 



Let 




v iat (t)=U (t)*v(t)U (t). 




Therefore, U v (t,s) = Uo(t)Ui nt (t,s)Uo(s)* satisfies the integral equation 



U v (t,t )=U (t-t )-i 



U (t-s)v(s)U v (s,to)ds 



(66) 



Denoting by y° the density operator of the crystal at time t — 0, the dynamics of 
the system is governed by the evolution equation 
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Y(t)=U v (t,0)fU v (t,Of. (67) 

Note that the conditions y° £ 5f (L 2 (R 3 )) and < y° < 1 are automatically propa- 
gated by d67l ). 

Considering v(f ) as a perturbation of the time-independent Hamiltonian H® el , and 
y(f ) as a perturbation of the ground state density operator y^ el ., it is natural to follow 
the same strategy as in the time-independent setting and introduce 

Using d66l ), ([67}, and the fact that y£ a is a steady state of the system in the absence 
of perturbation (Uo(t )"$ er Uo{t )* — }p er ), an easy calculation shows that Q(t) satisfies 
the integral equation 

Q(t) = Uo(t)Q(0)U {ty-iJ Uo(t-s)[v(s),^ + Q{s)}U (t -s)*ds. (68) 

We now assume that y° = y^ er , i.e. Q(0) — 0, and write (formally for the moment) 
Q(t ) as the series expansion 

G(0 = £(M0, (69) 

n=l 

where the operators Q n . v (t) are obtained, as in the time-independent case, by iden- 
tifying terms involving n occurrences of the potential v. In particular, the linear 
response is given by 

Qi,v(0 - -if U (t-s) [v(i),i£J U (t-s)*ds, (70) 
and the following recursion relation holds true 

Vn>2, Qn, v (t) = -i[ U (t-s)[v(s),Q n ^ hv (s)}U (t-s)*ds. (71) 

JO 

It is proved in [ 1 1 that for any n > 1 and any t > 0, the operator 2„,,.(f ) in d69l > 
belongs to and satisfies 

Vv/eL 2 (R 3 ), (v/|e n ,v(0lr>L2=0. 

In particular, Ti'o(Q n . v (t)) = 0- Besides, there exists b € R + such that for all f > 

ft ft\ ft,,- 1 

\\Qn,v{t)U<b n / .../ \\p{h)\\ L 2^...\\p{t n )\\ L 2^dt n ...dh, 
JO JO JQ 
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and there exists T > such that the series expansion d69l l converges in Js? uniformly 
on any compact subset of [0, T). Lastly, T = +°° if p e L°°(R + ,L 2 (R 3 ) n ¥). 

As in the time-independent setting, the frequency-dependent dielectric properties 
of the crystal can be obtained from the linear response dTOb . by defining the time- 
dependent independent-particle polarization operator 



Xo : L 1 (R,v c (L 2 (R 3 )n*'))->L"(R,L 2 (R 3 )n« : 

v PQt.v 



(72) 



and the time-dependent operators Jz? = — %qV c , e — v c (l + j£f)v c , £ = v c (l + 
— 1 v^ - 1 , and e = v c 1,/2 £Vc^ 2 . Due to the invariance of the linear response with 
respect to translation in time, all these operators are convolutions in time. In addi- 
tion they are ^-periodic in space. They can therefore be represented by frequency- 
dependent Bloch matrices [r K K /(a),q)], with K, K' in , q £ r* and 03 £ M. 
The Adler- Wiser formula states that the (electronic contribution of the) frequency- 
dependent macroscopic dielectric permittivity is given by the formula 

VkeR 3 , k T re M ((o)k = lim ( n |k f — — Y 

The mathematical study of this formula and of its possible derivation from rigorous 
homogenization arguments, is work in progress. 

We finally consider the self-consistent Hartree dynamics defined by 



Q(t) = U (t)Q°U (ty -ij U (t- s) [v(s) + v c (p QW ), & + Q(sj\ U (t - s)*ds, 

(73) 

for an initial condition Q° £ J(T, and for an external potential v(t) = v c {m(t)), where 
m(t) £ L 2 (R 3 ) ntf for all t. The solution Q(t) of is such that 7(f) = y£ er + Q(t) 
satisfies, formally, the time-dependent Hartree equation 



{Pylt) 



^per 



■m(t))*\-\-\y(t) 



The following result [ 10] shows the well-posedness of the nonlinear Hartree dy- 
namics. 



Theorem 3. Let m £ C l (R + ,L 2 (R 3 )n^). Then, for any Q° £ X, the time-dependent 
Hartree equation ( 1731 l has a unique solution in C°(R+,=S). Besides, for all t > 0, 
Q(t)£jrandTr (Q(t))=Tr (Q°). 
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Appendix: trace-class and self-adjoint operators 

It is well-known that any compact self-adjoint operator A on a separable Hilbert 
space can be diagonalized in an orthonormal basis set: 

A = £^|<fc)<&|, (74) 
1=1 

where (<j>i\<j>j) = dtj, and where the sequence (Ai),>i of the (real) eigenvalues of A, 
counted with their multiplicities, converges to zero. We have formulated d74l using 
again Dirac's bra-ket notation. The conventional mathematical formulation for d74l i 
reads 

1=1 

A compact self-adjoint operator A is called trace-class if 

i=i 

The trace of A is then defined as 

Tr(A):=£A, = £( e ,|Ah), 
i=i i=i 

the right-hand side being independent of the choice of the orthonormal basis (<?, ),> i . 
Note that if A is a non-negative self-adjoint operator, the sum Y,tTi( e i\A\ e d makes 
sense in M + U {+°°} and its values is independent of the choice of the orthonormal 
basis (c/)i>i- We can therefore give a sense to Tr(A) for any non-negative self- 
adjoint operator A, and this number is finite if and only if A is trace-class. 

The notion of trace-class operators can be extended to non-self-adjoint operators 
|j3Tll33ll . but we do not need to consider this generalization here. 

By definition, a compact operator A is Hilbert-Schmidt if A* A is trace-class. A 
compact self-adjoint operator A on Jff decomposed according to (f74l i is Hilbert- 
Schmidt if and only if 

I>| 2 <~. 

i>l 

Obviously any trace-class self-adjoint operator is Hilbert-Schmidt, but the converse 
is not true. 

In this contribution, we respectively denote by @i and &2 the spaces of trace- 
class and Hilbert-Schmidt operators acting onL 2 (R 3 ). We also denote by ^(L 2 (IR 3 )) 
the vector space of the bounded self-adjoint operators on L 2 (M 3 ). 

A classical result states that if A is a Hilbert-Schmidt operator on L 2 (M 3 ), then 
it is an integral operator with kernel in L 2 (R 3 x R 3 ). This means that there exists a 
unique function in L 2 (R 3 x R 3 ), also denoted by A for convenience, such that 



The Microscopic Origin of the Macroscopic Dielectric Permittivity of Crystals 37 
V</>eL 2 (R 3 ), (A0)(r)= / A(r,r')(j)(r')dr'. (75) 

Conversely, if A is an operator on L 2 (R 3 ) for which there exists a function A e 
L 2 (R 3 x R 3 ) such that (75$ holds, then A is Hilbert-Schmidt. 

If A is a self-adjoint Hilbert-Schmidt operator on L 2 (R 3 ) decomposed according 
to (1741 1. then its kernel is given by 

A(r/) = £A,<Mr)<Mr')- 

;>i 

If, in addition A is trace-class, then the density p^, defined as 

P/ ,(r) = X>#(r)| 2 , 

is a function of L 1 (R 3 ) and it holds 

Tr(A) = VXi= / p A (r)dr. 

For convenience, we use the abuse of notation which consists in writing Pa (r) = 
A(r,r) even when the kernel of A is not continuous on the diagonal {r = r'} C R 6 . 



References 



1. Adler, S.L.: Quantum theory of the dielectric constant in real solids. Phys. Rev. 126, 413^120 
(1962) 

2. Ambrosio, L., Friesecke, G. Giannoulis, J.: Passage from quantum to classical molecular 
dynamics in the presence of Coulomb interactions. Commun. Part. Diff. Eq. 35, 1490-1515 
(2010) 

3. Ambrosio, L., Figalli, A., Friesecke, G. Giannoulis, J., Paul, T.: Semiclassical limit of quan- 
tum dynamics with rough potentials and well posedness of transport equations with measure 
initial data. larXiv: 100 6.5388 yl (2010) 

4. Ananfharaman, A., Cances, E: Existence of minimizers for Kohn-Sham models in quantum 
chemistry. Ann. I. H. Poincare-An 26, 2425-2455 (2009) 

5. Arnold, A.: Self-consistent relaxation-time models in quantum mechanics, Commun. Part. 
Diff. Eq., 21(3-4), 473-506 (1996) 

6. Baroni, S., Resta, R.: Ab initio calculation of the macroscopic dielectric constant in silicon. 
Phys. Rev. B 33, 7017-7021 (1986) 

7. Cances, E, Deleurence, A., Lewin, M.: A new approach to the modelling of local defects in 
crystals: the reduced Hartree-Fock case. Commun. Math. Phys. 281, 129-177 (2008) 

8. Cances, E, Deleurence, A., Lewin, M.: Non-perturbative embedding of local defects in crys- 
talline materials. J. Phys.: Condens. Mat. 20, 294213 (2008) 

9. Cances, E, Lewin, M.: The dielectric permittivity of crystals in the reduced Hartree-Fock 
approximation. Arch. Ration. Mech. Anal. 197, 139-177 (2010) 

10. Cances, E, Stoltz, G.: in preparation 

1 1 . Catto, I, Le Bris, C, Lions, P.-L. : On the thermodynamic limit for Hartree-Fock type models. 
Ann. I. H. Poincare-An 18, 687-760 (2001) 



38 



Eric Cances, Mathieu Lewin and Gabriel Stoltz 



12. Dautray, R. and Lions, J.-L. Mathematical Analysis and Numerical Methods for Science and 
Technology. Vol. 5. Evolution Problems I. Springer- Verlag Berlin (1992) 

13. Dreizler, R., Gross, E.K.U.: Density functional theory. Springer Verlag, Berlin (1990) 

14. Engel, G. E., Farid, B.: Calculation of the dielectric properties of semiconductors. Phys. Rev. 
B 46, 15812-15827 (1992) 

15. Frank, R.L., Lieb, E.H., Seiringer, R., Siedentop, H.: Mlillers exchange-correlation energy in 
density-matrix-functional theory. Phys. Rev. A 76, 052517 (2007) 

16. Gajdos, M., Hummer, K., Kresse, G., Furthmliller, J., Bechstedt, E: Linear optical properties 
in the projector-augmented wave methodology. Phys. Rev. B 73, 0451 12 (2006) 

17. Gravejat, P., Lewin, M, Sere, E.: Ground state and charge renormalization in a nonlinear 
model of relativistic atoms. Commun. Math. Phys. 286, 179-215 (2009) 

18. Hainzl, C., Lewin, M., Sere, E.: Existence of a stable polarized vacuum in the Bogoliubov- 
Dirac-Fock approximation. Commun. Math. Phys. 257, 515-562 (2005) 

19. Hainzl, C, Lewin, M., Sere, E., Solovej, J. P.: A minimization method for relativistic electrons 
in a mean-field approximation of quantum electrodynamics. Phys. Rev. A 76, 052104 (2007) 

20. Hainzl, C, Lewin, M., Solovej, J.P.: The mean-field approximation in Quantum Electrody- 
namics: the no-photon case. Commun. Pur. Appl. Math. 60(4), 546-596 (2007) 

21. Hohenberg, P., Kohn, W.: Inhomogeneous electron gas. Phys. Rev. 136, B864-B871 (1964) 

22. Hybertsen, M.S., Louie, S.G.: Ab initio static dielectric matrices from the density-functional 
approach. I. Formulation and application to semiconductors and insulators. Phys. Rev. B 35, 
5585-5601 (1987) 

23. Hybertsen, M.S., Louie, S.G.:Ab initio static dielectric matrices from the density-functional 
approach. II. Calculation of the screening response in diamond, Si, Ge, and LiCl. Phys. Rev. 
B 35, 5602-5610 (1987) 

24. Kohn, W., Sham L. I.: Self-consistent equations including exchange and correlation effects. 
Phys. Rev. 140, A1133 (1965) 

25. Kunc, K., Tosatti, E.: Direct evaluation of the inverse dielectric matrix in semiconductors. 
Phys. Rev. B 29, 7045-7047 (1984) 

26. Lieb E.H.: Variational principle for many-fermion systems. Phys. Rev. Lett. 46, 457^159 
(1981) 

27. Lieb, E.H.: Density Functional for Coulomb systems. Int. J. Quantum Chem. 24, 243-277 
(1983) 

28. Pazy A.: Semigroups of Linear Operators and Applications to Partial Differential Equations, 
vol. 44 of Applied Mathematical Sciences. Springer, New York (1983) 

29. Pick R. M., Cohen, M.H., Martin R. M.: Microscopic theory of force constants in the adiabatic 
approximation. Phys. Rev. B 1, 910-920 (1970) 

30. Reed, M., Simon, B.: Methods of Modern Mathematical Physics. II. Fourier Analysis, Self- 
Adjointness. Academic Press, New York (1975) 

31. Reed, M., Simon, B.: Methods of Modern Mathematical Physics. IV. Analysis of Operators. 
Academic Press, New York (1978) 

32. Resta, R., Baldereschi, A.: Dielectric matrices and local fields in polar semiconductors. Phys. 
Rev. B 23, 6615-6624 (1981) 

33. Simon, B.: Trace ideals and their applications, vol. 35 of London Mathematical Society Lec- 
ture Note Series. Cambridge University Press, Cambridge (1979) 

34. Solovej, J.P.: Proof of the ionization conjecture in a reduced Hartree-Fock model. Invent. 
Math. 104, 291-311 (1991) 

35. Thomas, L.E.: Time dependent approach to scattering from impurities in a crystal. Commun. 
Math. Phys. 33, 335-343 (1973) 

36. Wiser, N.: Dielectric constant with local field effects included. Phys. Rev. 129, 62-69 (1963) 

37. Zhislin, G.M., Sigalov, A.G.: The spectrum of the energy operator for atoms with fixed nuclei 
on subspaces corresponding to irreducible representations of the group of permutations. Izv. 
Akad. Nauk SSSR Ser. Mat. 29, 835-860 (1965) 



